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Abstract. "The Price of Robustness" by Bertsimas and Sim [4j repre- 
sented a breakthrough in the development of a tractable robust counter- 
part of Linear Programming Problems. However, the central modeling 
assumption that the deviation band of each uncertain parameter is sin- 
gle may be too limitative in practice: experience indeed suggests that 
the deviations distribute also internally to the single band, so that get- 
ting a higher resolution by partitioning the band into multiple sub-bands 
seems advisable. The critical aim of our work is to close the knowledge 
gap about the adoption of a multi-band uncertainty set in Robust Opti- 
mization: a general definition and intensive theoretical study of a multi- 
band model are actually still missing. Our new developments have been 
also strongly inspired and encouraged by our industrial partners, which 
have been interested in getting a better modeling of arbitrary distribu- 
tions, built on historical data of the uncertainty affecting the considered 
real-world problems. In this paper, we study the robust counterpart of a 
Linear Programming Problem with uncertain coefficient matrix, when a 
multi-band uncertainty set is considered. We first show that the robust 
counterpart corresponds to a compact LP formulation. Then we inves- 
tigate the problem of separating cuts imposing robustness and we show 
that the separation can be efficiently operated by solving a min-cost flow 
problem. Finally, we test the performance of our new approach to Ro- 
bust Optimization on realistic instances of a Wireless Network Design 
Problem subject to uncertainty. 
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1 Introduction 

A fundamental assumption in classical optimization is that all data are exact. 
However, many real- world problems involve data that are uncertain or not known 
with precision, because of erroneous measurements or adoptions of approximated 
numerical representations. If such uncertainty is neglected, optimal solutions 
computed for nominal data values may become costly or infeasible. As a con- 
sequence, considering the impact of uncertainty on an optimization model is a 
critical issue when dealing with real-world problems. 

During the last years, Robust Optimization (RO) has become a valid method- 
ology to deal with optimization problems subject to uncertainty. A key concept 
of RO is to model uncertainty as hard constraints, that are added to the orig- 
inal formulation of the problem. This restricts the set of feasible solutions to 
robust solutions, i.e. solutions that are protected from deviations of the data. 
Such a robust approach is crucial when dealing with high risk events, such as 
aircraft scheduling |19j . or sensor placement in contaminant warning systems 
for water distribution networks [23]. In such settings, standard approaches like 
deterministic optimization or Stochastic Programming fail to protect against 
severe deviations, leading to unpredictable consequences. For an exhaustive in- 
troduction to the theory and applications of RO, we refer the reader to the book 
by Ben-Tal et al. [2], to the recent survey by Bertsimas et al. [3] and to the Ph.D. 
Thesis |8]. 

An approach to model uncertain data that has attracted a lot of attention 
is the so called -T-scenario set, introduced by Betsimas and Sim (BS) [1] and 
then adapted to several applications. The uncertainty model for a Linear Pro- 
gram (LP) considered in BS assumes that, for each coefficient a we are given 
a nominal value a and a maximum deviation d and that the actual value lies 
in the interval [a — d,a + d]. Moreover, a parameter r is introduced to repre- 
sent the maximum number of coefficients that deviate from their nominal value. 
Hence, r controls the conservativeness of the robust model and its introduction 
comes from the natural observation that it is unlikely that all coefficients devi- 
ate from their nominal value at the same time. A central result presented in BS 
is that, under the previous characterization of the uncertainty set, the robust 
counterpart of an LP corresponds to a linear formulation. This counterpart has 
the desirable properties of being purely linear and, above all, compact, i.e. the 
number of variables and constraints is polynomial in the size of the input of the 
deterministic problem. 

The use of a single deviation band may greatly limit the power of modeling 
uncertainty. This is particularly evident when the probability of deviation sensi- 
bly varies within the band: in this case, neglecting the inner-band behaviour and 
just considering the extreme values like in BS may lead to a rough estimate of the 
deviations and thus to unrealistic uncertainty set, which either overestimate or 
underestimate the overall deviation. Having a higher modeling resolution would 
therefore be very desirable. This can be accomplished by breaking the single 
band into multiple and narrower bands, each with its own r. Such model is 
particularly attractive when historical data about the deviations are available, 
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a very common case in real-world problems. A multi-band uncertainty set can 
indeed effectively approximate the shape of the distribution of deviations built 
on past observations, guaranteeing a much higher modeling power than BS. This 
observation was first captured by Bicnstock and taken into account to develop 
an RO framework for the special case of Portfolio Optimization [5] . It was then 
extended to Wireless Network Design |6ll0j . Yet, no definition and intensive the- 
oretical study of a more general multi-band model applicable in other contexts 
have been done. The main goal of this paper is to close such gap. 

We remark that investigating the adoption of a multi-band uncertainty set 
and studying the theoretical properties of the resulting model have been dis- 
cussed with and strongly encouraged by our industrial partners, such as British 
Telecom Italia (BT) and Nokia Siemens Networks (NSN), in past and present 
collaborations about real network design |6I7I10I15] . NSN, in particular, has been 
particularly interested in finding refined models for taking into account the ar- 
bitrary and non-symmetrical distributions of traffic uncertainty characterizing 
nation-wide optical networks |23j . Our modeling and theoretical developments 
are thus also strongly based on practical industrial needs. A better modeling of 
the traffic uncertainty affecting telecommunication networks has been also one 
of the critical objectives of our research activities in the German industrial re- 
search project ROBUKOM, aimed at developing new models and algorithms for 
the design of robust and survivable networks [7115] . in collaboration with NSN 
[2"T] and the German National Research and Education Network (DFN) [T2] . 

Contributions and Outline. In this work, we study the robust counterpart 
of an LP with uncertain coefficient matrix, when a multi-band uncertainty set is 
considered. The main original contributions are: 

1. a compact formulation for the robust counterpart of an LP; 

2. an efficient method for the separation of robustness cuts (i.e., cuts that 
impose robustness), based on solving a min-cost flow instance; 

3. computational experiments comparing the performance of solving the com- 
pact formulation versus a cutting plane approach on realistic wireless net- 
work design instances. 

In Section [5J we show that the robust counterpart of an LP under multi-band 
uncertainty corresponds to a compact Linear Programming formulation. We then 
proceed to study the separation problem of robustness cuts in Section [31 Finally, 
in Section [4j we test the performance of our new model and solution methods 
to Robust Optimization, to tackle the uncertainty affecting signal propagation 
in a set of realistic DVB-T instances of a wireless network design problem. 

1.1 Model and Notation 

In our study about the robust counterpart of Linear Programming Problems 
subject to multi-band uncertainty, we refer to a generic Linear Program of the 
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form: 

max Cj xj (LPP) 

tty < b{ i £ / 

> j£ J 



where / = {1,...,™} and J = {1, ...,n} denote the set of constraint and 
variable indices, respectively. 

We assume that the value of each coefficient a,j is uncertain and is equal 
to the summation of a nominal value ciij and a deviation lying in the range 
[dfj , <i^ + ], where dftj , £ K represent the maximum negative and positive 
deviations from a^, respectively. The actual value aij thus lies in the interval 
[aij + dfj , + dfj + ]. We note that assuming that the uncertainty only affects 
the elements of coefficient matrix does not limit the generality of our study, as 
uncertainty on the cost c and on the r.h.s. b can be included in a very straight- 
forward way in the coefficient matrix [3J. 

We derive a generalization of the Bcrtsimas-Sim uncertainty model by parti- 
tioning the single deviation band \dfj~ ,dfj + ] of each coefficient a,j into K bands, 
defined on the basis of K deviation values: 

-co < df; <■■• < eT 2 < d^ 1 < £a = < dh < cL 2 , < • • ■ < df A < +oo. 

Through these deviation values, we define: 1) a set of positive deviation bands, 
such that each band k £ {1, . . . , K + } corresponds to the range (d^f 1 , d^j]; 2) a 
set of negative deviation bands, such that each band k £ {K~ + 1, . . . , — 1, 0} 
corresponds to the range (d^f 1 , d^} and band k = K~ corresponds to the single 

value dfj (the interval of each band but k = K~ is thus open on the left). 
With a slight abuse of notation, in what follows we indicate a generic deviation 
band through the index k, with k £ K = {K~ , ... ,—1,0,1, ... , K + } and the 
corresponding range by (c^j -1 , d^]. 

Additionally, for each band k £ K we define a lower bound Ik and an upper 
bound Uk on the number of deviations that may fall in k, with Ik, Uk £ Z satisfy- 
ing < Ik < Uk < n. In the case of band k = 0, we assume that uq = n, i.e. we do 
not limit the number of coefficients that take their nominal value. Furthermore, 
we assume that X^fcgK < n , so that there always exists a feasible realization 
of the coefficient matrix. 

We remark that, in order to avoid an overload of the notation, we assume that 
the number of bands K and the bounds lk,Uk are the same for each constraint 
i £ I. Anyway, it is straightforward to modify all presented results to take into 
account different values of those parameters for each constraint. 

All the elements introduced above define what we call a multi-band uncer- 
tainty set Sm- We now proceed to study the robust counterpart of (LPP) under 
multi-band uncertainty. 
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2 A Compact Robust LP Counterpart 

The robust counterpart of (LPP) under a multi-band uncertainty set defined by 
Sm can be equivalently written as: 

max Cj Xj 
jeJ 

a%j Xj + DEVi(x, Sm) < h i E I 

jeJ 

xj > j G J, 

where DEVi(x, Sm) is the maximum overall deviation allowed by the multi-band 
uncertainty set for a feasible solution x when constraint i is considered. Note 
that we replace the actual value of a coefficient a%j with the summation of the 
nominal value a%j and a deviation falling in exactly one of the K bands. The 
computation of DEVi(x,Sm) corresponds to the optimal value of the following 
pure 0-1 Linear Program (note that in this case the index i is fixed): 

DEVi(x, S M ) = max E E d% Xj y% (DEVOl) 

j£j k£K 

h < y% ^ u k k e K (1) 

3&J 

E 4 < 1 3 e J ( 2 ) 

k£K 

y-j^ {0,1} jeJ,keK. (3) 

The binary variables y\j indicate if the deviation of a coefficient lies in band 
fc. The constraints © ensure that each coefficient deviates in at most one band 
(actually these should be equality constraints, but, for assumption uq = n made 
in Scction ll.il we can consider inequalities). Finally, the constraints (TTJ) impose 
the upper and lower bounds on the number of deviations falling in each band fc. 
Thus, an optimal solution of (DEVOl) defines a distribution of the coefficients 
among the bands that maximizes the deviation w.r.t. the nominal values, while 
respecting the bounds on the number of deviations of each band. 

We now prove that the robust counterpart of (LPP) under a multi-band uncer- 
tainty set Sm can be reformulated as a compact Linear Program. To this end, 
consider first the linear relaxation of (DEVOl): 



max V" ^ 



d% xj y% (DEVOl-RELAX) 

jeJ fee A' 

i k < y% - Uk ke K ( 4 ) 

E y% ^ 1 3 e j (5) 

keK 

Vij > jeJ,ke K, (6) 



6 C. Busing, F. D'Andreagiovanni 



where we dropped constraints y^ < 1 since they are dominated by constraints 
©• 

Proposition 1. The polytope described by the constraints of (DEV01-RELAX) 
is integral. 

Proof. We start by rewriting all the constraints of (DEV01-RELAX) into the 
form a T y < j3 obtaining the following matrix form: 



A yi 



-I 


-/ 




\ 

-/ 




I 


I 




I 






1 ■ ■ • 1 




1...J 





Vn 



\yfn J 



( : \ 

-Ik 



< 



Ilk 



\ ■■ J 



Consider now the sub-matrix Di obtained from Di by eliminating the top layer 
of blocks (— 7| — I\ ■ ■ ■ | — 7). It is easy to verify that Di is the incidence matrix 
of a bipartite graph: the elements of the two disjoint set of nodes of the graph 
are in correspondence with the rows of the two distinct layers of blocks in Di. 
Moreover, every column has exactly two elements that are not equal to zero, one 
in the upper layer and one in the lower layer. Being the incidence matrix of a 
bipartite graph, Di is a totally unimodular matrix |20j . 

In order to show that also the original matrix Di is totally unimodular, we 
first need to recall the equivalence of the following three statements [20]: 1) A 
is a totally unimodular matrix; 2) a matrix obtained by duplicating rows of A 
is totally unimodular; 3) a matrix obtained by multiplying a row of A by -1 is 
totally unimodular. Since Di can be obtained from Di by duplicating each row 
of the upper block, and multiplying each row of the duplicated block by -1, Di 
is totally unimodular. 

As Di is totally unimodular and the vector gi is integral, it is well-known 
that the polytope defined by Diyi < gi and yi > is integral, thus completing 
the proof. □ 



Thanks to this integrality property, we can exploit strong duality to prove the 
main result of this section. 
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Theorem 1. The robust counterpart of (LPP) under a multi-band uncertainty 
set Sm is equivalent to the following compact Linear Program: 

max Cj Xj (RLP) 

Y x i - Y lk v i + Y uk w i + Y bi ieI 

je.J keK keK jeJ 

-v* + Wi + z\ > dtj xj i g l,j e J,k G K 

v%, w\ > iel,ke K 

z{ > i G I,j G J 

> J G J. 

Proof. As first step, consider the dual problem of (DEV01-RELAX): 



min ~h + Y Uk w i + z i (DEV01-RELAX-DUAL) 

k£K k£K j£.J 

-v\ + + z{ > d% Xj j G J,k G K 



wf, < > fee Jv 



i e J, 



where the dual variables , , z\ are respectively associated with the primal 
constraints fll © of (DEV01-RELAX) defined for constraint i. 

Since (DEVO 1-REL AX) is feasible and bounded by definition of the scenario 
uncertainty set Sm, also (DEV01-RELAX-DUAL) is feasible and bounded and 
the optimal values of the two problems are the same (strong duality). Then, 
by Proposition 1, we can replace the maximum total deviation DEVi(x,SM) 
with problem (DEV01-RELAX-DUAL), obtaining the compact Linear Program 
(RLP). This concludes the proof. □ 

In comparison to (LPP), this compact formulation uses 2 • K • m + n • m addi- 
tional variables and includes K ■ n ■ m additional constraints. Similar reasonings 
can be also done to derive a compact robust counterpart of a Mixed-Integer 
Linear Program with uncertain coefficient matrix. 



3 Separation of Robustness Cuts 

In this section, we consider the problem of testing whether a solution x G K n 
is robust feasible, i.e. J2j£j®v x j + D EVi(x, Sm) < bi for every scenario S G 
Sm and i G I. This problems becomes important if we adopt a cutting plane 
approach instead of directly solving the compact robust counterpart (RLP): we 
start by solving the nominal problem (LPP) and then we check if the optimal 
solution is robust. If not, we generate a cut that imposes robustness (robustness 
cut) and we add it to the problem. This initial step can then be iterated as in a 
typical cutting plane method [2D] . 
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In the case of the Bertsimas-Sim model, the problem of separating a robust- 
ness cut is extremely simple [13]: given a solution x G ffi™, for each constraint 
i 6 J, the problem consists of sorting the deviations d^ + Xj in non-increasing 
order and choose the highest P, deviations. If for some i the sum of these de- 
viations exceeds bi — cLaXj then we have found a robustness cut to add. 
Otherwise, x is feasible and robust. 

In the case of multi-band uncertainty, this simple approach does not guaran- 
tee the robustness of a computed solution. However, we prove that for a given 
solution x G W 1 and a constraint i <E I, checking the robust feasibility of x cor- 
responds to solving a min-cost flow problem pQ and thus can be done efficiently. 
The min-cost flow instance associated with the robustness check is denoted by 
(G, c) x and defined as follows. G is a directed graph that contains one vertex 
Vj for each variable index j G J, one vertex Wk for each band k G K and 
two vertices s, t that are the source and the sink of the flow. So the set of ver- 
tices is V = {Jj eJ {vj} U {J keK {wk} U {s,t}. The set of arcs A is the union of 
three sets A\, A%, A3. A\ contains one arc from s to every variable vertex Vj, i.e. 
A\ = {(s, Vj) : j G J}. A2 contains one arc from every variable vertex Vj to every 
band vertex Wk, i.e. A2 = {(vj, Wk) ■ j G J, k G A"}. Finally, A3 contains one arc 
from every band vertex Wk to the sink t, i.e. A3 = {(wk, t) : k G K}. By construc- 
tion, G(V, A) is bipartite and acyclic. Each arc a G A is associated to a triple 
(l a , u a ,c a ), where l a , u a are lower and upper bounds on the flow that can be sent 
on a and c a is the cost of sending one unit of flow on a. The values of the triples 
(l a , u a , c a ) are set in the following way: (0, 1, 0) when a £ A\\ (0, 1, — d^xj) when 
a = (vj, Wk) G A2] {Ik, Uk, 0) when a = (wk, t) G A3. Finally, the amount of flow 
that must be sent trough the network from s to t is equal to n. We denote by F$ 
the set of feasible integral flows of value n of the min-cost flow instance (G, c) l x . 
The cost of an (s, i)-flow / G Fi is defined by c(f) = J2 a eA c afa- We recall that 
an integral min-cost flow of value n can be computed in polynomial time, using 
for example the successive shortest path algorithm [T|. 

We now prove that by solving the min-cost flow instance (G, cf x defined 
above, we obtain the maximum deviation for a constraint i and a solution x. 

Theorem 2. Let x G R™ and letSu be a multi-band uncertainty set. Moreover, 
let (G, c) l x be the min-cost flow instance corresponding with x and a constraint 
i E I of (LPP) and built according to the previously presented rules. 
The solution x is robust feasible w.r.t. Sm for constraint i if and only if 

a -a; — c*(x) < 6, 

where c*(x) is the minimum cost of a flow of the instance (G, c) x . 

Proof. Before proceeding to the core of the proof, we need to make some obser- 
vations. Let Yi be the set of feasible solutions to problem (DEV01), i.e.: 

Yi = {y, G {0, 1} |J| |K| : Vl satisfies ©-©}. 

Each of these solutions represents a feasible assignment of the coefficients of con- 
straint i to the deviation bands K. Moreover, let d(x,yt) be the total deviation 
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associated with a vector x £ K™ and a feasible assignment yi GYi, i.e.: 

d(x, yi ) = E E d % x j y%- 

We note that this is actually the objective function of problem (DEV01). 

It is easy to verify that for each feasible assignment that distributes less than 
n = \J\ coefficients among the bands, there exists a feasible assignment that 
distributes all the n coefficients among the bands and entails at least the same 
total deviation, i.e.: 

v v\ e E E v$ < n 

d(x,yf)>d(x,yj). 

Thanks to the last observation, in order to find a feasible assignment y^ £ Yi 
that maximizes the total deviation d(x, yi), we can concentrate attention on the 
subsets % C Yi of assignments which distribute all the n coefficients among the 
bands: 

Yi = {Vi G {0, 1} 1 J| W ■■ y l satisfies ©-© and ^ ^ y% = n). 

jeJkeK 

We now prove the lemma by showing that, given a solution x £ K™ and a 
constraint i £ /, there exists a one-to-one correspondence between a feasible 
assignment yi € Yi with total deviation d{x,yi) and an integral flow / of value 
n and cost c(f) — —d(x,yi) of the min-cost flow instance (G,c) l x . 

Once that this correspondence is proved, we can solve the problem of finding 
an assignment with maximum deviation DEVi (x,Sm) by finding an integral flow 
of minimum cost of the instance (G, cf x . Indeed the following chain of equations 
holds: 

DEV^x^Sm) = max d(x, yi) = — min — d{x,yi) = — min c(f) = —c*{x). 

y.eY, yiGYi f£Fi 



Given a feasible assignment yi £ Yi with total deviation d(x,yi), we construct 
an integral flow / for the min-cost flow instance (G,c) x , in the following way: 

f{s.v 3 ) = 1 Vj £ J (7) 

f( Vj , Wk ) = 1 <=> Vij = 1 Vj £ J, k £ K (8) 

iW) = E4- VfceA ^- ( 9 ) 

It is clear that for each j £ J, k £ K the flow components /( S)tl .), f( Vj ,w k ) respect 
the capacity of the corresponding arcs. In the case of f{ Wkl t)i the capacity of arc 
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(wk,t) is respected as © holds and Ik < Ylj^jVij < Uk, Vfe G K by definition 
of yi. By (JTJ), the value of the flow / is YljeJ f(s.vj) = \ J\=n and its cost is: 

c(/) = ~4> ^ /(«*>™*) = ~ d % X J y% = - d ( x >Vi)- 

Concerning flow conservation, by definition of ?/,;, we know that each coefficient 
of the constraint i is assigned to exactly one band, i.e. 3 unique k G K : yfj = 1, 
Vj £ J. As a consequence, flow conservation is respected in every vertex Vj with 
j G J. In the case of a vertex u>fc with k £ K, we have: 

/ (tiifc.t) = 2/y = f(vj,w k )i 

jeJ jeJ 

so the flow conservation is respected. Finally, in the case of the sink t, we have 

an ingoing flow of J2keK f(w h ,t) = E je j J2keK Vij = n > which is thus cc l ual to 
the value of the flow /. 



Given a feasible integral flow / G Fi of value n and cost c(/) of the min-cost 
flow instance (G,c) l x , we construct a deviation assignment G {0, 1}' J ' ' A ' , in 
the following way: 

w&- = 1 /fe.^o = l Vj g J, fc g X. (10) 

Since the value of / is n by definition and since each arc incident to a vertex Vj , 
j G J has lower and upper capacity bounds of and 1, respectively, i) the flow 
on each arc (s, Vj) is unitary, ii) the ingoing flow of each vertex Vj is unitary, 
and by flow conservation we have: 

Y Y = Y Y y*j = n - 

j£j k£K jeJ k£K 

Moreover, by flow conservation and by the integrality of /, the outgoing flow of 
vj is unitary and is sent to a single node Wk- Therefore: 

Y hi,™*) = Y y % = 1 

keK keK 

Finally, by definition of lower and upper capacity bounds of an arc ft Wk ,t) an d 
by flow conservation, for every k G K we have: 

^fc < f(w k ,t) < u k, 

f(w k ,t) = Y] f (Vj,Wk) = ^ Vij 

jeJ je.J 

h < < u k 

j£j 
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By resuming all the properties of yt that we proved, it follows that y^ G Y± and 
the associated total deviation is equal to: 

-d(x,yi) = E E ~ d ij x i y% = E E ~ d ij x i hv 3 ,w k ) = c(/). 

This concludes the proof of the existence of a one-to-one correspondence between 
the set Yi of feasible deviation assignments and the set Fi of integral flows 
and of the relation DEVi(x,Sm) = ~ c i{ x ) between the optimal values of the 
maximum deviation problem (DEV01) and the min-cost flow problem of the 
instance (G, c) l x . The statement of the problem follows from this correspondence. 

□ 

From Theorem [21 we can immediately derive the following corollary: 

Corollary 1. A solution x £ W] is robust feasible w.r.t. Sm if o,nd only if 

a-x — c*{x) < bi Mi 6 /, 

where c*(x) is the minimum cost of a flow of the instance (G, c) l x . 

According to this corollary, we can test the robustness of a solution x € R™ by 
computing an integral flow f of minimum cost c*(x) in (G,c) l x for every i <E I 
(we recall that this can done in polynomial time). If a'x — c*(x) < bi for every 
i, then x is robust feasible. Otherwise x is not robust and there exists an index 
i E I such that a'x — c*(x) > bi and thus: 

je.J jeJ kcK 

is valid for the polytope of the robust solutions and cuts off the solution x. 

4 Computational Study 

In this section, we test our new modeling and solution approaches to Robust 
Optimization on a set of realistic instances of the Power Assignment Problem, 
a problem arising in the design of wireless networks. In particular, we compare 
the efficiency of solving directly the compact formulation (RLP) with that of a 
cutting plane method based on the robustness cuts presented in Section [3] In the 
case of the Bertsimas-Sim model, such comparison led to contrasting conclusions 
(e.g., [T3IT4]). 

The Power Assignment Problem. The Power Assignment Problem (PAP) 
is the problem of dimensioning the power emission of each transmitter in a 
wireless network, in order to provide service coverage to a number of user, while 
minimizing the overall power emission. The PAP is particularly important in 
the (re)optimization of networks that are updated to new generation digital 
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transmission technologies. For a detailed introduction to the PAP and the general 
problem of designing wireless networks, we refer the reader to |18ll0lllfT7j . 

A classical LP formulation for the PAP can be defined by introducing the fol- 
lowing elements: 1) a vector of non-negative continuous variables p that represent 
the power emissions of the transmitters; 2) a vector p max of upper bounds on p 
that represent technology constraints on the maximum power emissions; 3) a ma- 
trix A of the coefficients that represent signal attenuation (fading coefficients) for 
each transmitter-user couple; 4) a vector of r.h.s. 5 (signal-to-interference thresh- 
olds) that represent the minimum power values that guarantee service coverage. 
Under the objective of minimizing the overall power emission, the PAP can be 
written in the following matrix form: 

min l'p s.t. Ap > S, < p < P max (PAP) 

where exactly one constraint a\p > 5i is introduced for each user i to represent 
the corresponding service coverage condition. 

Each entry of matrix A is classically computed by a propagation model and 
takes into account many factors (e.g., distance between transmitter and receiver, 
terrain features). However, the exact propagation behavior of a signal cannot be 
evaluated and thus each fading coefficient is naturally subject to uncertainty. Ne- 
glecting such uncertainty may provide unpleasant surprises in the final coverage 
plan, where devices may turn out to be uncovered for bad deviations affecting 
the fading coefficients (this is particularly true in hard propagation scenarios, 
such as dense urban fabric). For a detailed presentation of the technical aspects 
of propagation, we refer the reader to [22] . 

Following the ITU recommendations (e.g., [IB]), we assume that the fading 
coefficients are mutually independent random variables and that each variable 
is log-normally distributed. The adoption of the Bertsimas-Sim model would 
provide only a rough representation of the deviations associated with such dis- 
tribution. We thus adopt the multi-band uncertainty model to obtain a more 
refined representation of the fading coefficient deviations. In what follows, we 
denote the Bertsimas-Sim and the multi-band uncertainty model by (BS) and 
(MB), respectively. 

Computational Results. In this computational study, we consider realistic 
instances corresponding to region-wide networks that implement the Terrestrial 
Digital Video Broadcasting technology (DVB-T) [TB] and were taken as reference 
for the design of the new Italian DVB-T national network. The uncertainty set is 
built taking into account the ITU recommendations [16] and discussions with our 
industrial partners in past projects about wireless network design. Specifically, 
we assume that each fading coefficient follows a log-normal distribution with 
mean provided by the propagation model and standard deviation equal to 5.5 
dB [16| . In our test-bed, the (MB) uncertainty set of a generic fading coefficient 
dij is constituted by 3 negative and 3 positive deviations bands (i.e., K = 6). 
Each band has a width equal to the 5% of the nominal fading value ay. Thus 
the maximum allowed deviation is +/- 0.15 • ay. For each constraint i, the 
bounds Ik , Uk on the number of deviations are defined considering the cumulative 
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distribution function of a log-normal random variable with standard deviation 
5.5 dB. The (BS) uncertainty set of each constraint considers the same maximum 
deviation of (MB) and the maximum number of deviating coefficients is r = 
[0.8 • u max ], where u max = max{u fe : k e A'\{0}}. This technically reasonable 
assumption on r ensures that (BS) does not dominate (MB) a priori. 

The computational results are reported in Table [TJ The tests were performed 
on a Windows machine with 1.80 GHz Intel Core 2 Duo processor and 2 GB 
RAM. All the formulations are implemented in C++ and solved by IBM ILOG 
Cplex 12.1, invoked by ILOG Concert Technology 2.9. We considered 15 in- 
stances of increasing size corresponding to realistic DVB-T networks. The first 
column of Table [1] indicates the ID of the instances. Columns |/|, |J| indicate the 
number of variables and constraints of the problem, corresponding to the number 
of user devices and transmitters of the network, respectively. We remark that 
the coefficient matrices tend to be sparse, as only a (small) fraction of the trans- 
mitters is able to reach a user device with its signals. Columns |/ + |, | J + \ indicate 
the number of additional variables and constraints needed in the compact robust 
counterpart (RLP). Columns PoR% report the Price of Robustness (PoR), i.e. 
the deterioration of the optimal value required to guarantee robustness. In par- 
ticular, we consider the percentage increase of the robust optimal value w.r.t. the 
optimal value of the nominal problem, in the multi-band case (PoR% (MB)) and 
in the Bertsimas-Sim case (PoR% (BS)). Column At% reports the percentage 
increase of the time required to compute the robust optimal solution under (MB) 
by using the cutting plane method presented in Section[3]w.r.t. the time needed 
to solve the compact formulation (RLP). Finally, column Protect% is a measure 
of the protection offered by the robust optimal solution and is computed in the 
following way: for each instance, we generate 1000 realizations of the uncertain 
coefficient matrix and we then compute the percentage of realizations in which 
the robust optimal solution is feasible. This is done for both (MB) and (BS). 

Looking at Table [IJ the first evident thing is that the dimension of the com- 
pact robust counterpart under (MB) is much larger than that of the nominal 
problem. However, this is not an issue for Cplex, as all instances are solved 
within one hour and in most of the cases the direct solution of (RLP) takes less 
time than the cutting plane approach (At% < 0). Anyway, for the instances 
of greater dimension the cutting plane approach becomes competitive and may 
even take less time (At% > 0). Concerning the PoR, we note that under (MB) 
imposing robustness leads to a sensible increase in the overall power emission, 
that is anyway lower than that of (BS) in all but two cases. On the other hand, 
such increase of (MB) is compensated by a very good 90% protection on average. 
In the case of the PAP, (MB) thus seems convenient to model the log-normal 
uncertainty of fading coefficients, guaranteeing good protection at a reasonable 
price. Moreover, though (BS) offers higher protection for most instances, it is 
interesting to note that the increase of Protect% of (BS) w.r.t (MB) is lower 
than the corresponding increase of PoR% of (BS) w.r.t (MB). 
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Table 1. Overview of the computational results 













-TOKyo 


foti/o 




rrotectyo 


rrotect/o 


ID 


l 1 


J 


IJ+I 


J + l 


(MB) 


(BS) 


At% 


(MB) 


(BS) 


Dl 


95 


153 


3519 


10098 


8.3 


10.1 


-18.7 


88.20 


92.53 


D2 


103 


197 


4728 


14184 


7.2 


9.4 


-19 


91.35 


92.47 


D3 


105 


322 


7406 


21252 


6.8 


8.8 


-16.9 


93.12 


96.40 


D4 


105 


473 


10406 


28380 


7.4 


7.2 


-15.1 


92.08 


91.42 


D5 


108 


569 


13087 


37554 


9.2 


11.4 


-13.6 


89.23 


90.29 


D6 


157 


1088 


27200 


84864 


6.6 


9.1 


-6.2 


85.46 


87.55 


D7 


165 


1203 


31278 


101052 


7.1 


9.5 


-4.9 


87.91 


89.16 


D8 


171 


1262 


32812 


106008 


8.7 


10.8 


-4.1 


89.40 


93.08 


D9 


178 


1375 


35750 


115500 


9.6 


10.2 


-2.8 


90.11 


91.90 


D10 


180 


1448 


39096 


130320 


7.9 


9.6 


-1.7 


91.54 


95.32 


Dll 


180 


1661 


46058 


159456 


7.2 


9.5 


0.6 


94.77 


96.70 


D12 


181 


1779 


49812 


170784 


7.5 


10.1 


1.8 


88.22 


90.16 


D13 


183 


1853 


53737 


189006 


8.1 


10.3 


3.3 


91.34 


92.21 


D14 


183 


1940 


56260 


197880 


10.3 


9.7 


3.1 


86.50 


85.18 


D15 


185 


2183 


63307 


222666 


8.4 


10.8 


4.1 


91.09 


92.70 



5 Conclusions and Future Work 

In this work, we presented new theoretical results about multi-band uncertainty 
in Robust Optimization. Surprisingly, this natural and practically very relevant 
extension of the classical single band model by Bertsimas and Sim has attracted 
very little attention and we have thus started to fill this theoretical gap. We 
showed that, under multi-band uncertainty, the robust counterpart of an LP is 
linear and compact and that the problem of separating a robustness cut can 
be formulated as a min-cost flow problem and thus be solved efficiently. Tests 
on realistic network design instances showed that our new approach performs 
very well, thus encouraging further investigations. Future research will focus on 
refining the cutting plane method and enlarging the computational experience 
to other relevant real-world problems. 
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